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INTRODUCTION 

Obsolete  or  unwanted  munitions,  rocket  propellants,  and  manufacturing  wastes  require 
treatment  at  Department  of  Defense  (DOD)  and  Department  of  Energy  (DOE) 
facilities.  One  of  the  most  widely-used  treatment  methods  is  open  burning  (OB)  and 
open  detonation  (OD)  of  the  material.  Currently,  the  material  destroyed  in  a  single 
detonation  generally  ranges  from  100  to  5000  lbs,  while  the  quantity  removed  by  a  burn 
can  be  larger  and  last  from  minutes  to  an  hour.  OB/OD  operations  are  confined  to 
daytime  with  atmospheric  stability  conditions  ranging  from  convective  or  highly  unstable 
to  near  neutral. 

OB /OD  activities  produce  air  pollutants  and  require  predictions  of  pollutant  concen¬ 
trations  to  obtain  an  operating  permit.  The  pollutants  include  SO2,  NOi,  particulates, 
volatile  organic  compounds  and  toxic  materials  such  as  metals,  semivolatile  organics, 
etc.*^  Large  detonations  also  generate  large  quantities  of  dust/soil  that  are  entrained  by 
the  rising  contaminant  cloud.  OB/OD  sources  differ  from  most  traditional  air  pollution 
sources  in  that  they  have:  1)  instantaneous  or  short-duration  releases  of  buoyant 
material  rather  than  continuous  releases,  and  2)  ambient  exposure  times  for  the  clouds 
that  can  be  much  less  than  the  typical  averaging  times  (  >  l.hr)  of  air  quality  standards. 

Atmospheric  dispersion  models  are  used  to  estimate  pollutant  concentrations  given 
information  on  the  source  and  meteorological  conditions.  However,  there  is  currently  no 
recommended  EPA  dispersion  model  to  address  the  unique  features  of  OB/OD  sources. 
The  most  commonly-used  approach  is  INPUFF, ^  a  Gaussian  puff  model,  but  this  has 
several  limitations  as  briefly  discussed  below  and  elsewhere.^  As  a  result,  a  model 
development  program  was  initiated  under  the  DOD/DOE  Strategic  Environmental 
Research  and  Development  Program. 

This  paper  briefly  summarizes  the  model  development  effort  which  is  divided  into 
“operational”  and  “research”  components.  In  the  following,  we  give  a  brief  discussion 
of  the  background  and  overall  model  design  and  then  describe  the  operational  model 
components  for  instantaneous  (OD)  and  short-duration  (OB)  sources;  the  OD  model  is  a 
Gaussian  puff  approach  whereas  the  OB  framework  consists  of  integrated-puff  and  plume 
models.  The  combined  OB/OD  model  includes:  1)  a  continuous  treatment  of  dispersion 
as  the  release  condition  varies  from  instantaneous  to  continuous,  2)  cloud  and  plume 
rise  obtained  from  appropriate  entrainment  models,  3)  cloud  and  plume  penetration 
of  elevated  inversions,  4)  relative  (puff)  and  total  dispersion  based  on  modern  scaling 
concepts  for  the  plane^ry  boundary  layer  (PBL),  and  5)  a  capability  for  the  use  of 
onsite  profiles  of  wind,  temperature,  and  turbulence  from  a  mobile  meteorological 
platform.  The  current  OB/OD  model  focuses  on  the  unstable  PBL. 
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BACKGROUND  AND  MODEL  DESIGN 


Background 

The  development  of  an  OB/OD  dispersion  model  has  considered:  1)  the  limitations  of 
existing  models.  2)  current  knowledge  of  turbulence  and  dispersion  in  the  PBL,  and  3)  a 
mobile  meteorological  platform  under  development. 

Limitations  of  existing  models.  INPUFF  has  been  used  to  model  OB/OD  sources 
and  can  handle  dispersion  from  individual  puffs  or  clouds  or  from  a  sequence  of  puffs  in 
a  short-duration  release.  Although  the  Gaussian  puff  approach  is  suitable  for  OB/OD 
sources,  INPUFF  has  the  following  limitations;  1)  It  adopts  dispersion  parameters 
(cTyiCTj)  from  the  Pasquill-Gifford  (PG)- curves.  2)  It  includes  Briggs’  (Ref.  4)  plume 
rise  expressions  which  apply  to  continuous  releases  rather  than  to  instantaneous  sources 
(puffs,  clouds)  and  does  not  address  thermal  penetration  of  elevated  inversions  capping 
the  PBL.  3)  It  assumes  Gaussian  velocity  statistics  for  the  turbulence,  whereas  the 
vertical  velocity  statistics  in  the  unstable  PBL  are  positively  skewed.®  The  skewness 
should  be  included  for  vertical  dispersion. 

For  OB/OD  sources,  the  PG  curves  are  deficient  in  that  they:  1)  are  bzised  on  dispersion 
from  a  ground-level  source  and  short  downwind  distances  (<  1  km),  and  2)  are  selected 
using  smface  meteorology,  which  does  not  account  for  the  PBL’s  vertical  structure.  For 
large  detonations,  source  buoyancy  can  carry  emissions  to  several  100  m  or  the  PBL  top; 
one  must  then  deal  with  dispersion  over  the  entire  PBL. 

PBL  turbulence.  Dispersion  in  the  PBL  depends  on  the  turbulence  length  and 
velocity  scales  which  differ  for  the  unstable  or  convective  boundary  layer  (CBL)  and 
the  stable  boundary  layer  (SBL).  For  the  CBL,  the  length  and  velocity  scales  are  the 
CBL  depth  h  and  the  convective  velocity  scale  iw*;  w*  =  {gwdohfQai  where  g  is  the 
gravitational  acceleration,  w6o  is  the  turbulent  heat  flux  at  the  surface,  and  0a  is  the 
ambient  potential  temperature.  Typical  values  of  w*  and  h  at  midday  over  land  are  1-2 
m/s  and  1-2  km.  Within  the  “mixed  layer  ’  (O.lh  <  z  <  h),  the  mean  wind  speed  and 
turbulence  components — longitudinal  cr^,  lateral  (7„,  and  vertical  cr-u, — vary  little  with 
height  z\  in  strong  convection,  (Xv^cTw  —  0.6u;.. 

For  the  SBL,  the  turbulence  is  much  weaker  with  eddy  sizes  proportional  to  z  near  the 
surface  and  typically  ~  10s  of  meters  or  less  in  the  upper  part  of  the  SBL.  Models  and 
observations  show  that  the  velocity  scale  is  the  friction  velocity  tt.  (Ref.  5),  which  is 
typically  ~  0.1  m/s  in  strong  stable  stratification. 

Knowledge  of  the  PBL  turbulence  structure  has  been  included  in  a  number  of  models  for 
air  quality  applications.®  One  example  is  AERMOD’^  for  industrial  source  complexes. 

Mobile  meteorological  platform.  A  mobile  meteorological  platform  is  being 
developed  at  NOAA’s  Environmental  Technology  Laboratory  to  obtain  the  PBL 
variables  necessary  for  modeling  since  many  DOD  facilities  are  in  remote  locations. 

The  platform  design  includes;  1)  a  radar  wind  profiler  for  obtaining  the  three  wind 
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components  up  to  ~  3  km.  2)  a  radio  acoustic  sounding  system  (RASS)  for  temperature 
measurements.  3)  a  mini-SODAR  for  measuring  winds  and  cTu,  to  a  height  of  ~  200  m. 

4)  a  mini-lidar  system  for  obtaining  the  PBL  depth  /i.  and  5)  a  porta.ble  meteorological 
station  for  measuring  near-surface  winds,  temperature,  turbulence,  and  heat  flux.  The 
dispersion  model  is  being  designed  for  efficient  use  of  these  measurements. 

Overall  Model  Design 

The  overall  model  design  includes:  1)  a  simple  computational  or  operational  framework 
for  routine  problems,  and  2)  a  more  detailed  or  research  model  for  nonroutine  problems. 
In  the  operational  approach,  a  Gaussian  puff  model  is  adopted  for  instantaneous  sources 
and  puff,  integrated-puff,  and  plume  models  for  short-duration  releases.  For  the  research 
framework,  a  Lagrangian  particle  and/or  puff  approach  is  planned.  Both  frameworks 
will  be  considered  for  “onsite”  use  in  a  real-time  operational  mode  using  data  from  the 
mobile  meteorological  platform,  i.e.,  for  day-to-day  decisions  on  OB/OD  operations.  The 
operational  puff  and  plume  models  would  be  used  for  climatological  analyses  needed  in 
risk  assessments. 

In  modeling,  the  important  aspects  to  address  are:  1)  all  source-related  features 
including  the  instantaneous  or  short-duration  nature  of  the  release,  buoyancy-induced 
rise  and  dispersion,  and  cloud  or  plume  penetration  of  eleyated  inversions,  2)  relative 
and  absolute  dispersion  expressions  that  explicitly  include  PBL  turbulence  variables,  3) 
meteorological  variables  including  their  vertical  profiles  from  the  mobile  platform,  and  4) 
a  treatment  for  puff  and  plume  dispersion  about  complex  terrain. 

The  models  discussed  in  the  following  sections  address  points  1  and  2  above  and  must 
be  expanded  to  include  points  3  and  4.  Further  development  also  will  address:  1)  a 
more  complete  description  of  initial  source  effects  (detonation  cloud  size  and  height) 
and  inversion  penetration,  2)  a  more  complete  PBL  turbulence  parameterization,  3) 
averaging  time  effects  on  concentration,  4)  the  entrained  dust  source  term,  and  5) 
deposition  of  gases  and  particles. 

INSTANTANEOUS  SOURCES 


Dispersion  Model  /nr  n  \ 

In  many  OB/OD  applications,  estimates  of  peak  ground-level  concentrations  (GLCs) 

are  required  for  averaging  times  ranging  from  seconds  or  minutes  to  an  hour  or  longer. 
In  making  such  estimates,  we  must  account  for  the  stochastic  or  random  nature  of 
turbulence  and  dispersion  in  the  PBL.  That  is,  we  must  recognize  that  the  observed 
concentration  at  a  receptor  is  a  random  variable  and  should  be  predicted  statistically 
through  a  probability  distribution.®’^  The  distribution  can  be  parameterized  by  a 
functional  form  such  as  a  gamma  or  clipped-normal  distribution®  and  requires  a 
minimum  of  two  variables  to  characterize  it— the  mean  concentration  C  and  the  root- 
mean-square  (rms)  concentration  fluctuation  Cc,  which  is  a  measure  of  the  width  of 
the  distribution.  The  peak  concentration  then  can  be  defined  by  a  specified  percentile 
value  of  the  cumulative  probability  distribution,  e.g.,  the  99.9th  percentile  level.  In  the 
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following,  we  discuss  approaches  for  predicting  C.  Gc-  cloud  rise,  and  the  dispersion 
parameters:  the  functional  form  of  the  distribution  remains  to  be  selected. 


Mean  concentration.  For  instantaneous  sources  or  detonations,  a  Gaussian  puff  model 
is  adopted  for  predicting  the  short-term  mean  concentration  (C)  field  relative  to  the  puff 
centroid.  The  C  is  the  expected  or  average  concentration  that  would  be  observed  if  the 
same  experiment — same  source  and  meteorological  conditions — were  repeated  a  large 
number  of  times.®  The  C  is  given  by: 


C  = 


Q 


(27r)^^^(T  ry^r 


exp 


{x-Utf 


{z  -  KY 

2^r"x 


(1) 


where  Q  is  the  pollutant  mass  released,  U  is  the  mean  v/ind  speed,  t  is  the  travel  time. 
he  is  the  effective  puff  height,  and  <7rx,  o’ryj  and  <7rz  are  the  puff  standard  deviations 
or  relative  dispersion  parameters  in  the  x,  y,  and  z  directions,  respectively.  Here,  he  = 
hs  +Ah{x)  where  hg  is  the  source  height  which  is  generally  zero  for  OB/OD  sources,  and 
Ah  is  the  cloud  rise  due  to  buoyancy;  x  and  y  are  the  distances  in  the  mean  wind  and 
‘crosswind  directions. 


The  maximum  mean  concentration  Cc  at  the  puff  centroid  is  given  by  Cc  = 
Q/[(27r)®/^crrx<J’ryO’rz]>  where  the  relative  dispersion  parameters  are  generally  different 
in  the  three  directions.  In  the  following,  we  assume  =  o'ry  =  <^rz  =  CTr- 


The  C  field  including  puff  meandering  is  given  by  Eq.  (1),  but  with  O' rx  ^  Ory^  Orz 
replaced  by  the  absolute  dispersion  parameters — o'x,o’y,<^z-  A  Gaussian  distribution 
for  C  is  applicable  in  the  SBL  where  the  probability  density  function  (p.d.f.),  p^,  of 
the  vertical  velocity  w  is  Gaussian.  However,  for  the  CBL,  a  skewed  w  p.d.f.  is  more 
consistent  with  laboratory  and  field  data.  A  skewed  p.d.f.  is  adopted  here  and  is 
parameterized  by  the  superposition  of  two  Gaussian  distributions.^® 

The  C  field  due  to  an  ensemble  of  meandering  puffs  is  derived  from  following  the 
same  approach  cis  applied  to  continuous  plumes.^®  The  resulting  expression  for  C  is 


Q 


exp 


(x  -  f/t)2 
2<^1 


{Z- he-  Zj)^ 


(2) 


where  ffzj  =  cr^jxlU  and  =  WjxJU  with  j  =  1,2.  The  Xj,  w],  and  cr«,j  {j  =  1,2) 
are  the  weight,  mean  velocity,  and  standard  deviation  of  each  Gaussian  p.d.f.  comprising 
Pt0.  Equation  (2)  applies  for  short  distances  such  that  the  plume_interaction  with  the 
ground  or  elevated  inversion  is  weak.  The  complete  expression  for  C  includes  multiple 
cloud  reflections  at  the  ground  and  PBL  top. 


The  time-averaged  concentration  can  be  found  from  the  dose  where  the  partial  dose  is 
defined  by  t/>(x,  y,  z,  t)  =  C(x,  y,  z,  t')dt'  and  the  total  dose  by  i/>oo  =  ^  ^)- 

For  clouds  with  short  passage  times  over  a  receptor,  the  average  concentration  C  can  be 
obtained  from  C  =  (t/»(<2)  -  ^{ti))/Ta,  wh^e  the  averaging  time  Ta  =  ta  -  If  the  puff 
passage  time  Aarx/U  is  less  than  Ta,  then  C  =  ^oc/^a- 
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Rms  concentration  fluctuation.  For  elevated  sources  such  as  a  buoyant  cloud,  the 
rms  concentration  fluctuation  at  the  ground  is  largest  close  to  the  source  and  decreases 
with  increasing  distance.^’^^  As  a  first  step  in  predicting  ctc,  we  adopt  Gifford’s^^ 
meandering  plume  or  cloud  model  which  applies  to  the  near-source  region;  specifically, 
it  applies  when  the  large-scale  eddies  in  the  PEL  cause  the  wandering  or  meandering 
of  a  small  growing  plume  or  cloud.  The  cloud  only  occasionally  reaches  the  ground, 
but  it  does  so  in  high  concentration  due  to  the  small  local  spread  dr'-,  thus,  it  creates 
large  “spikes"  in  the  GLC  distribution  and  a  large  concentration  variance.  In  the 
following,  we  model  the  mean  square  concentration  {c^),  where  the  angle  brackets  denote 
an  ensemble  average,  and  then  find  dc  from  dc  =  ((c^)  -  (see  Ref.  9). 


In  Gifford’s  model,  the  instantaneous  concentration  distribution  in  a  cloud  is  assumed 
to  be  a  nonrandom,  axisymmetric  Gaussian  distribution  about  the  instantaneous  cloud 

centroid  Xc  =  {xciVci  ^c)’- 


c(x,xc) 


Q 

(27r)^/^cr^ 


exp 


(x  -  Xc)^  (y  -  _  jz  -  Zc)'^ 

2d}  2d}  2d} 


(3) 


Concentration  fluctuations  arise  due  to  the  randomness  in  Xc,  Vc  and  which  is  caused 
by  the  meandering  of  the  cloud. 


The  (c^)  is  obtained  by  averaging  the  quantity  (P‘  from  Eq:  (3)  over  the  joint 
displacement  p.d.f.,  p{xci  0i  which  is  a  function  of  time.  In  the  well-mixed  CBL,  the 
turbulent  velocity  fluctuations  in  the  x,  y,  z  directions  are  uncorrelated  throughout  the 
bulk  of  the  layer;  likewise,  the  centroid  displacements  in  the  three  directions  should  be 
uncorrelated  random  variables.  Thus,  assuming  independence  of  the  displacements  in 
the  three  directions,  we  have  p(xc,  t)  =  Px{^c,i)  •  PyiVc^t)  -  Pz{zc^  0-  Using  Gifford  s 
approach,  the  (c^)  is  found  from 


C^{x;Xc)Px(Xc,  t)  ■  PyiVc,  t) 


'Vc 


Pz(^Zf2i  t')(iX(^ciyc(iZc  . 


(4) 


The  displacement  p.d.f.s  in  (4)  are  found  from  the  velocity  p.d.f.s — Pu(u),  Pu(u),  and 
Pw{w),  where  u,v,w  are  the  random  velocity  components  in  the  x,y,z  directions.  Close 
to  the  source,  we  assume  that  Xc  =  ut,  Pc  =  Zc  =  +  wt.  The  PxiPy  ^.nd  Pz  can 

then  be  derived  as  in  Weil.^®  Taking  Pu  and  p„  as  Gaussian,  one  finds 


Px(^c) 


-exp 


(Xc  -  {Xc))^ 


2al 


■  ,  PyiVc)  = 


-exp 


vl ' 


(5) 


where  (Xx  =  (xc)  =  Ut,  and  dy  =  dyt.  For  p^,,  a  bi-Gaussian  p.d.f.  or  superposition 
of  two  Gaussian  distributions  is  adopted,  which  results  in  a  similar  distribution  for  p*  as 


noted  earlier  (Eq.  2). 


By  substituting  the  derived  Px,  Py,  and  Pz  into  Eq.  (4)  and  carrying  out  the  integration, 
we  obtain 


(c^)  = 


{2T:)^d}dxe(^ye 


exp 


( 


(x  —  Ut)'^ 
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where  Oxe.  —  (<^r  +  (^ye  =  (cTr  "h  and  cr^ej  =  (cr^  +  with  j  =  1  or 

2.  The  Aj.  zj,  and  are  the  same  variables  as  in  Eq.  (2).  The  above  {(?)  expression 
applies  to  short  distances  such  that  the  cloud  interaction  with  the  ground  and  elevated 
inversion  is  weak;  similar  to  Eq.  (2),  we  include  reflections  at  z  =  0,/i  to  obtain  the 
complete  expression  for  (c^). 

Note  that  Eq.  (6)  is  of  the  same  form  as  Eq.  (2)  except  for  the  omission  of  the  2  in  the 
denominator  of  the  exponential  terms.  This  is  consistent  with  the  {<?)  expression  for 
passive  releases  in  Gaussian  turbulence  (see  Refs.  11  and  12). 

Cloud  rise  and  inversion  penetration.  Scorer^^  combined  theory  and  laboratory 
experiments  to  obtain  the  following  expression  for  cloud  rise  in  a  neutral  environment 

Ah  =  2.35(MTt  +  •  (7) 


Mr  and  Ft  are  the  initial  momentum  and  buoyancy  of  the  cloud  and  are  given  by 


Att  o 

Mr  =  -Trriwo  and  Ft  = 


qQt 

CpPa^a 


(8) 


where  Wo,  To-,  and  Qt  are  the  initial  velocity,  radius,  and  heat  content  of  the  thermal, 

Cp  is  the  specific  heat  of  air,  and  Pa  and  ©a  are  the  ambient  air  density  and  potential 
temperature.  The  Qt  can  be  determined  from  the  mass  of  the  detonation  and  its  heat 
content,  H  =  1100  kcal/kg  TNT  equivalent. 

Scorer  also  found  the  puff  radius  to  be  r  =  aA/it,  where  Aht  is  the  cloud  top  height  and 
a  is  an  empirical  entrainment  coefficient,  q  ranged  from  0.14  to  0.5  with  a  mean  of  0.25. 
The  relative  dispersion  ar  —  rf  \/2. 

Using  field  observations,  Weil^^  confirmed  that  Eq.  (7)  was  a  good  fit  to  data  over  a 
wide  range  of  times.  Thus,  Eq.  (7)  is  suitable  for  the  rise  of  a  cloud  before  it  is  limited 
by  stable  stratification,  e.g.,  an  elevated  inversion. 

For  cloud  penetration  of  an  elevated  density  jump,  results  have  been  found  from 
laboratory  experiments  in  a  nonturbulent  environment.  Richards obtained  an 
empirical  expression  for  the  fraction  P  of  the  cloud  penetrating  the  jump:  P  ~  1  - 
Q.5Api/ApTi,  where  Ap<  is  the  density  jump  and  Apn  is  the  average  density  excess  of 
the  cloud  when  it  reawihes  the  jump.  This  can  be  applied  to  atmospheric  detonations 
and  elevated  temperature  inversions  by  relating  Ap  to  the  cloud^temperatiire  exce^  A©, 
expressing  A©  in  terms  of  Ft  and  r,  and  evaluating  the  resulting  expression  at  z  —  h. 
From  Ap/ Pa  =  AG/&a  and  the  above  approach,  one  finds  A©  =  {3/Air)QT/{PaCp‘^^)- 
Richards’  expression  can  then  be  rewritten  as  P  =  1  —  {2w/Z){AQiPaCpa^h^ /Qt)^  where 
A©i  is  the  temperature  jump.  The  hg  is  assumed  to  be  zero  so  that  the  cloud  radius  at 
z  =  his  ah.  The  above  relationship  shows  the  strong  sensitivity  of  P  to  ah. 

A  more  realistic  temperature  distribution  above  the  CBL  is  a  constant  dQa/^^- 
Experiments  simulating  this  distribution  as  well  as  a  jump  above  a  well-mixed  layer  are 
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currently  underway  in  a  salt-stratified  tank  at  the  EPA  Fluid  Modeling  Facility  in  North 
Carolina. 


Dispersion  parameters.  For  clouds,  cr^  is  dominated  by  entrainment  for  short  times 
with  <7r  =  (Trb  =  O.I8A/1.  wherc  subscript  b  denotes  buoyancy-induced  spread.  At 
intermediate  times  {t  <  Ti),  the  cr^  may  be  dominated  by  ambient  turbulence  in  the 
inertial  subrange  with  <Tr  ~  cTra  =  where  T^,  is  the  Lagrangian  time^scale, 

e  is  the  turbulent  kinetic  energy  dissipation  rate,  oi  is  a  constant  (see  Thomson  ),  and 
subscript  a  denotes  dispersion  due  to  ambient  turbulence.  At  long  times  (t  »  T^),  cr^a  = 
{2(rljTLty^^  for  homogeneous  isotropic  turbulence.  For  ffra^  we  use  an  interpolation 
expression  of  the  form  Cra  —  +  02^/Tc,)  to  satisfy  the  intermediate-  and 

long-time  results.  In  addition,  e  can  be  written  as  c  =  hcr^/Ti,  in  homogeneous  isotropic 
turbulence. 

In  a  strong  CBL,  the  following  approximations  can  be  made  for  z  >0.lh:  c  ~  OAwl/h, 
Cyj  0.6u;.,  and  Tl  ~  0.7hlw»  (Ref.  10).  These  approximations  coupled  with  e  = 
bal/TL  lead  to  6  =  0.78.  To  satisfy  the  long-time  ara  limit,  we  must  have  ca  =  0.62ai; 
Oi  is  estimated  to  be  0.57  from  Thomson’s  two-particle  model  results.  The  resulting 
parameterization  for  ara  in  the  CBL  is 


^ra 

~h 


0.36X3/2 

1-F0.51X 


with 


tw.x 

""  Hh 


(9) 


where  we  have  assumed  t  =  x/U  and  X  is  a  dimensionless  distance  or  travel  time. 


To  connect  the  short-,  intermediate-,  and  long-time  relative  dispersion  regimes  in  a 
continuous  manner,  we  adopt  the  following  parameterization;  a^  =  a%  4-  a^a-  For  clouds 

dominated  by  buoyancy,  arh  —  QA2Frp  t^/2. 

The  total  or  absolute  dispersion  is  necessary  to  estimate  the  C  for  a  meandering  puff  or 
plume.  The  Cx  and  ay  in  Eq.  (2)  can  be  obtained  from  a  parameterization  of  Taylor's 
theory:  ax  =  o-ut/(l  +  tl2TLx)^^'^  and  similarly  for  ay.  The  Tlx  is  the  Lagrangian  time 
scale  for  the  u  component  and  can  be  parameterized  by  T^x  ^  (^u/h^  ate.  (e.g.,  see  Ref. 
6).  For  the  CBL  and  the  results  below,  we  use  T^x  =  T^y  =  0.7/i/u;*  and  cr„  =  cr„  = 

0.6tt;,. 


Example  Results 

We  have  computed  the  Cc  in  the  buoyant  puflF  and  the  mean  GLC  along  y  =  0  due 
to  a  meandering  puff  for  0.1  W  <  50  tons.  The  <Tr,  o'*,  and  Cy  were  calculated  as 
described  in  the  previous  section.  In  the  following,  the  cloud  buoyancy  is  characterized 
by  its  dimensionless  value 

in  the  examples  below,  we  use  =  2  m/s,  h  =  1000  m.  and  f/  =  5  m/s. 


Figure  1  shows  calculated  values  of  the  peak  (Cc)  SO2  concentrations  in  a  detonation^ 
cloud.  The  cloud  SO2  mass  is  estimated  from  Q  —  W  •  Ef,  where  Ef  (=  2.23  x  10  ; 
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Andrulis^)  is  the  SO2  emission  factor.  At  small  x.  all  of  the  curves  have  the  same  slope; 

Cc  oc  because  Crj,  oc  and  the  Cr  is  dominated  by  <7^5  near  the  source.  For 

200  m  <  X  <  2000  m,  some  curves  exhibit  a  short  region  of  a  nearly  constant  Cc 
with  x;  this  is  due  to  puff  trapping  in  the  CBL  and  a  temporary  limitation  on  vertical 
dispersion  due  to  the  elevated  inversion.  At  large  distances  (x  >  10  km),  clouds  for  all 
cases  become  uniformly  mixed  in  the  vertical  but  continue  to  spread  laterally.  Thus,  the 
Cc  tends  towards  Q/cr^^  oc  Q/x\  the  curves  for  W  —  10-50  tons  show  that  Cc  varies 
inversely  with  x  for  large  distances. 

Figure  2  shows  the  mean  GLC  C  along  the  puflf  centerline  (y  =  0)  for  the  same  range 
of  W  and  Ft  values  as  in  Fig.  1.  This  mean  is  for  an  ensemble  of  meandering  puffs  and 
is  obtained  from  Eq.  (2)  with  reflection  terms  included.  Several  interesting  features  are 
found;  1)  A  non- monotonic  variation  occurs  in  the  maximum  GLC  Cm  with  W  and  Ft- 
2)  The  variation  in  Cm  for  0.1  <  W  <  50  tons  is  only  about  a  factor  of  4  even  though 
the  range  in  Q  is  a  factor  of  500;  the  weak  dependence  on  Q  is  attributed  to  the  increase 
in  A/i  with  Ft-  3)  The  Cm  is  of  the  order  of  0.1  yg/m^,  which  is  the  lower  bound  for  Cc 
in  Fig.  1.  4)  The  increase  in  the  distance  to  the  maximum  concentration  with  W  is  due 
to  the  increase  in  A/i  with  Ft  or  W. 

For  (Tc  >  C,  the  mean  GLC  along  y  =  0  due  to  an  ensemble  of  meandering  puffs 
probably  has  little  to  do  with  the  observed  centerline  GLC  in  an  individual  puff.  This  is 
attributed  to  the  large  variability  in  individual  concentration  observations  when  cTc  >  C. 
As  noted  earlier,  the  computed  C  in  Fig.  2  would  be  used  together  with  a  modeled  <Tc  in 
a  concentration  probability  distribution  to  estimate  the  peak  short-term  GLC  that  could 
occur  downstream  of  the  detonation. 

Some  preliminary  calculations  of  ctc/C  at  the  distance  of  the  maximum  GLC  were  made 
for  the  examples  in  Fig.  2.  For  W  =  0.1, 1,  and  5  tons,  the  corresponding  adC  was  6,  2, 
and  1.8;  for  W  =  10  to  50  tons,  the  aJC  <  1  and  was  decreasing  rapidly  with  increasing 
distance.  The  rapid  decrease  was  attributed  to  the  large  cloud  size  (cTr)  relative  to  h  and 
the  approach  of  the  concentration  distribution  to  a  vertically  well-mixed  value.  The  Oc 
model  requires  generalization  to  address  the  Cc  at  large  distances;  x  >  (2  -  Z)Uhfw^, 
or  about  5  -  7.5  km  in  the  current  example.  This  will  be  pursued  in  the  future. 

SHORT-DURATION  RELEASES 


Dispersion  Model 

For  short-duration  releases  5or  bums,  our  general  approach  is  an  integrated  puff  model  in 
which  the  short-term  mean  concentration  relative  to  the  puff  centerline  is 


Qr  dt' 

(27r  cx^vy^Tz 


(11a) 


/  =  exp 


(x  -  U{t  -  t'))^  _  _  {z  -  he)^' 

20-2^  2(t2^  2(t2^  . 


(116) 
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where  x  =  (x.j/.z).  is  the  puff  emission  time,  tr-is  the  total  release  duration.  Qr  is 
the  continuous  source  emission  rate,  cTn  =  o’rxii  ~  t')^  similarly  for  CFry,crrz-^  The 
integration  in  (11a)  can  be  carried  out  analytically  for  limiting  forms  of  C7rx{t  —  t),  etc., 
but  must  be  done  numerically  in  general. 

The  integrated  puff  model  also  is  used  for  estimating  the  mean  concentration  due  to  a 
sequence  of  meandering  puffs  by  replacing  the  relative  dispersion  (cTn,  etc  in  Eq.  11)  by 
the  absolute  dispersion — cTi,  ay,  a^. 

In  the  following,  we  calculate  the  Cc  for  a  short-duration  release  by  numerically 
integrating  Eq.  (11)  with  an  —  =  ar  =  CL\€}^^{t  —  /(IT  ~  ^  )ITl)- 

Similarly,  we  compute  the  GLC  C  for  the  same  release  but  for  the  meandering  puffs  with 
=  (Ty  =  0.6w.(t-f')/(l+0-5(i-<')/TLx)^/^  and  Tlz  =  O.T/i/w*.  For  an  infinitely-long 
release,  the  Cc  and  C  values  should  reduce  to  those  for  a  continuous  plume  as  shown 
below;  thus  the  plume  results  serve  as  an  upper  bound  to  the  concentration  values  for 
the  short-duration  release. 


The  mean  concentration  field  relative  to  a  plume  centerline  is  given  by 

^  Qr  (  _  {z-he?\ 

^  -  2-KUaryan^'^^\  2a^y  2a^,  )  ' 


(12) 


Here,  the  plume  rise  is  attributed  to  buoyancy  and  is  given  by  Ah  —  jU 

and  its  radius  is  r  =  0.4 Ah  (Ref.  17).  Source  momentum  effects  can  be  included  in  the 
future.  As  with  the  puff  model,  we  will  assume  ary  =  arz  =  (^r  and  a^  =  a^^,  +  cr,.^. 

The  ara  is  given  by  Eq.  (11)  and  the  plume  arb  =  r/\/2  =  /U .  The  Cc  = 

Qr/{2irUa^)  from  Eq.  (12);  these  expressions  are  expanded  to  include  reflection  at  z  = 

0,h. 

Example  Results 

Results  are  presented  for  the  dimensionless  concentration  CcUh^ /Qr  for  the  plume  and 
integrated-puff  models,  with  reflection  at  z  =  0,  h  included  in  both.  The  buoyancy  of  the 
continuous  source  is  characterized  by  the  dimensionless  buoyancy  flux  (Ref.  10). 


F,  = 


Fb 

Uwlh  ' 


(13) 


where  Fb  is  the  source  buoyancy  flux  which  is  proportional  to  the  source  heat  flux. 

Figure  3  shows  the  CcUh^/Qr  as  a  function  of  X  for  the  plume  model  with  F»  in  the 
range  0.001  <  F,  <  0.3.  The  large  variation  in  the  dimensionless  Cc  at  short  range  {X  < 
1)  is  due  to  the  buoyancy-induced  dispersion  <7^5.  As  can  be  seen.  CcUh^/Qr  decreases 
systematically  and  significantly  with  an  increase  in  F,  due  to  the  increase  in  arb  with 
F..  For  A  >  1.  the  curves  converge  to  the  same  limit  because  at  long  times  the  ar  is 
dominated  by  a^,  which  is  independent  of  F.. 


Figure  4  presents  the  CcUh'^/Qr  for  both  the  plume  and  integrated-puff  models  for  F.  = 
0.01  and  various  values  of  tr.  =  trW./h.  the  dimensionless  release  duration.  The  time 
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scale  h/vj,  =  500  s  for  the  w.  (=  2  m/s)  and  h  (=  1000  m)  used  here,  so  that  tr  ranges 
from  50  s  to  500  s  or  about  1  to  8  min.  As  can  be  seen,  the  plume  result  (solid  curve)  is 
an  upper  bound  to  the  integrated-puff  model  results.  The  X  value  corresponding  to  the 
departure  of  the  integrated-puff  solution  from  the  plume  solution  increases  as  tr»  does. 

For  X  >  5,  the  integrated-puff  concentrations  can  be  significantly  less  than  the  plume 
concentrations;  this  is  due  to  the  finite  duration  of  the  source. 

Figure  5  shows  the  dimensionless  mean  GLC,  CUh^/Qr,  for  the  meandering  plume 
and  for  a  sequence  of  meandering  puffs  (finite-duration  release).  The  trends  with  X 
and  tr*  are  similar  to  those  in  Fig.  4,  but  the  near-source  variation  of  CUh^fQr  with 
X  differs  from  that  for  Cc  due  to  the  diflferent  dispersion  rates  in  relative  (Fig.  4)  and 
absolute  dispersion  (Fig.  5).  Note  that  at  A"  =  0.1,  the  GLC  value  is  about  an  order 
of  magnitude  smaller  than  the  Cc  (Fig.  4)  because  the  ground  is  significantly  removed 
from  the  elevated  plume  centroid.  However,  at  X  =  10,  the  concentration  values  are 
quite  similar;  this  occurs  because  the  concentration  distribution  becomes  uniform  in  the 
vertical  due  to  plume  trapping  and  the  horizontal  spreads  (relative  and  absolute)  are 
about  the  same. 
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Figure  1.  SO2  concentration  at  detonation  cloud  centroid  as  a  function  of  downwind 
distance  and  detonation  mass  W. 


Figure  2.  Mean  ground-level  SO2  concentration  along  puflf  centerline  as  a  function  of 
downwind  distance  and  detonation  mass;  see  Fig.  1  for  key  to  lines. 
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Figure  3.  Dimensionless  concentration  at  plume  centroid  as  a  function  of  dimensionless 
downwind  distance  and  dimensionless  buoyancy  flux  F.. 


NOTE  TO  EDITORS 

Under  the  new  federal  copyright  law, 
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Fisure  4.  Diineiisioiiless  concentration  at  plume  or  puff  centroid  as  a  function  of 
dimensionless  downwind  distance  and  dimensionless  release  duration 


X 

Figure  5.  Dimensionless  mean  ground-level  concentration  along  plume  or  puff  centerline 
as  a  function  of  dimensionless  downwind  distance  and  dimensionless  release  duration 
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